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Factorial moments and cumulants are usually defined with respect to the unconditioned 
Poisson process. Conditioning a sample by selecting events of a given overall multiplicity 
N necessarily introduces correlations. By means of Edgeworth expansions, we derive 
^ ■ generalized cumulants which define correlations with respect to an arbitrary process rather 

than just the Poisson case. The results are applied to correlation measurements at fixed 
N, to redefining short-range vs. long-range correlations and to normalization issues. 
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In high energy multiparticle production, the variation of correlations with overall multiplicity 
N is most sensitive to the underlying dynamics |Lj, |], |3], f|, || . This sensitivity provides a good 
opportunity to test and discriminate between various models. However, with the notable ex- 
ceptions of Refs. H [7|], little attention has so far been devoted to a systematic treatment of 
iV-dependent correlation measures. 

Selecting events by fixed values of N is an example of conditioning. One of the merits of this type 
of conditioning is the separation of so-called "short-range" from "long-range" correlations. The 
latter arise from non-poissonian overall multiplicity distributions as well as the iV-dependence 
of the conditional single particle density p(x\N) (where x is any phase space variable of interest). 

Every conditioning, however, introduces new correlations which can be regarded as "unphys- 
ical" since they are merely a consequence of the selection procedure. Consequently, the very 
concept of a correlation must be rethought and generalized for the case of conditioned samples. 

In this paper, we seek to clarify and generalize factorial moments and cumulants for the case 
of fixed-multiplicity samples. This generalization suggests a natural strategy for separating 
nontrivial fixed-iV correlations from overall multiplicity conditioning effects. Generalized cu- 
mulants will hence also be useful in the study of pion interferometry as a function of N and 
other situations requiring conditioning. 

Multiparticle final states are best understood in terms of point processes [§, ||. The concept 
of correlation in point processes is usually based on the fundamental Poisson process, which is 
characterized by the factorial moment generating functional (f.m.g.f 1.) 



where p(x) = {1/ ai)da\ nc \/ dx is the differential inclusive cross-section at point x. Equivalently, 
for the Poisson process every finite-dimensional counting distribution, i.e. the joint probability 
of the particle counts n = . . . , Um) in an arbitrary finite partition of an overall phase space 



where 7(n m ; u m ) = exp(—i/ m )i/^ / n m \ with mean multiplicity v m = J Qm p(x)dx. 

Particles created in a Poisson process are completely uncorrelated in the sense that all factorial 
cumulant densities of order q > 1 vanish, 



while Ki(x) = p(x). 

By contrast, the multinomial process is defined by placing a fixed number iV of particles indepen- 
dently of each other into the domain Q tot according to a probability density p(x\N) = p(x\N)/N\ 






m=l 
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it has a f.m.g.f 1. 
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and its finite-dimensional counting distributions are multinomials. 



Despite the fact that the particles are placed independently of each other, the standard facto- 
rial cumulants of the multinomial process are nevertheless nonzero, thereby indicating purely 
"external" correlations: 



.mult , 



x 1: ...,x q \N) = (q-l)\ (-N) 1 -* J] p(xi\N) . 
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In other words, standard cumulants measure not only the "internal" correlations between 
particles but also the deviation of the overall multiplicity N from the Poisson distribution 
7 (jV; N) = exp(-N)N N /N\. 

It thus seems natural to ask for quantities that are sensitive not to such "external" deviations 
from the overall Poissonian but only to "internal" correlations, i.e. deviations from the multi- 
nomial. 



The answer, we believe, lies in the use of generalized cumulants which successfully separate the 
two. To justify their existence and form, we turn to the alternative definition of cumulants as 
coefficient functions in an Edgeworth expansion. In most of the literature, the latter is under- 
stood as the expansion of a continuous univariate distribution P(x) around the gaussian (cf. 
IU|). For integer counts, it is more natural to define a multivariate discrete Edgeworth expan- 



sion around the Poisson process. In its most general form, the discrete Edgeworth expansion 
of a point process / with respect to the Poisson process 7 reads, in terms of their f.m.g.f l.'s 
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with x = (xi, . . . ,x q ) and X(x) = X(xi) ■ ■ ■ X(x q ). The univariate expansion of the overall 
multiplicity distribution f(N) in terms of the Poissonian is 



f(N) = exp 



00 ttrf 

:4-A0(-v) + E^r(-v)< 

g=2 I' 



(7) 



where V is the discrete difference operator such that V-y(n; v) = 7(71; v) — -y(n — I; v). 



The extension to a multivariate expansion for the joint counting distributions in arbitrary finite 
partitions of Q tot then reads 



f(n) = exp 



M 



Wi (m) - v m ] (-V m ) + , 
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where V m 7(n; V) = j(n; v) — H m > r y(n m '—5 m f m ; v), and fi^(m) is the g-fold integral of K q (x) over 
the domains Q mi , . . . , Q mq . 

It can be verified easily that the coefficients in these expansions are identical with the gth 
order multivariate factorial cumulants of the process /. Therefore standard cumulants can be 
defined as the coefficient functions of the Edgeworth expansion around the Poisson process. 

This definition can be generalized in an obvious way. From Eq. (|6]), it is easy to show that the 
generating functional Qf of a process / can be expanded in terms of the f.m.g.fl. of any other 
process hft 

Q f [X) = exp £ L_2. / dx[nf(x) - n h q (x)]X(x) 

q=l 1- J 

where again n{ and represent the standard factorial cumulants of / and h. 



Qh[X], (9) 



It is quite natural, then, to define generalized cumulants as the coefficients in the expansion of 
/ in terms of h, 

«//* = «/-«} (10) 

for q > 2, while for q=l, we set by convention If h is the Poisson process, these 

generalized cumulants reduce to the standard cumulants. 

The generalized cumulants can be found as functional derivatives of the ratio of generating 
functionals, 



i=i 



6X(xi)J \ Q h 



(11) 

A=0 

where the Poisson generating functional Q 1 has to satisfy (5<5 7 /5A)|a=o = K i- This ensures 
that the above convention for K^ h is satisfiedP]. Note that K^ h is generally a standard cumulant 
to some well defined process only if the process / is divisible by h in the sense that the ratio 
Qf/Qh itself again satisfies all requirements for a f.m.g.fl. 

The general validity of Eq. (||) means that generalized cumulants will work for any type of con- 
ditioning, i.e. show up the real correlations of a process / with respect to a given conditioned 
reference process h. 

We now apply the concept of generalized cumulants to fixed-iV multiplicity samples. The 
discussion will range from the definition of internal cumulants for fixed- N measurements to 
resulting cumulant formulae for fixed-bin and correlation integral measurements, followed by 
a short discussion of improved short-range vs. long-range correlation formulae, summing over 
windows in N, internal moments, and normalization. 

Internal cumulants are the generalized cumulants with respect to the multinomial process (§): 

njjSlN) = K q {x\N) - k^\x\N) , (12) 

2 It is understood that such relations are valid only for processes having well-defined cumulants of all orders. 

3 One can view the ratio Qh/Q-y as the central moment generating functional of the process h. 
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with K™ ult (x\ N) given by Eq. (|5|) and K q (x\N) the standard cumulant measured for a fixed- N 
process. These internal cumulants will vanish everywhere, even on a differential level, if the 
N particles in the full domain Q tot are placed independently of each other following the same 
probability density p(x\N), i.e. behave multinomially. 

Univariate internal cumulants are obtained from the multivariate ones by integrating all q 
variables Xi,...,x q over the same domain Q m , Kg(m|iV) = J Qm dx Kg(x\N), yielding 

4(m\N) = (n^) N -3{n^} N (n m ) N + 2(l-N- 2 ){n m } 3 N , 

etc., with tv® = n\/(n—q)\. Event averages ()jv are taken over subsamples of fixed N. These 
integrated internal cumulants differ from the inclusive ones [|TTJ by a correction factor (1 — iV 1-9 ) 
in the last term. 

Integrated internal cumulants of order q > 2 vanish for any fixed-N distribution if the inte- 
gration domain Q m is enlarged to the full domain Qtot, in contrast to the standard cumulants 
which integrate to (— l) 9-1 (q — 1)! N. 

The same procedure can be used in the correlation integral method. In the notation of Ref. 
0, if a = J2j ©(e — Xjj) is the "sphere count" around particle i in event a and (b)^ = 
(%2j 0(e — X°-j)) N is the average of counts within other events b, the lowest-order internal 
cumulants are correspondingly given by 



4(m = (E^-a-^ 1 )^]} , as) 

\ i In 

4(e\N) = /j2[a [2] -(b^) N -2a(b) N + 2(l-N-')(b)l}), (14) 
\ t In 

and so on; the same (1 — N l ~ q ) factors appear as coefficients of the last terms {b)^ 1 for higher- 
order internal cumulants. Again, these internal correlation integrals vanish for distances e large 
enough to encompass the full domain f2 to t- 

Clearly, the multinomial correction will be most important for small iV and lowest order. The 
correction may or may not be substantial, depending on the numerical size of the uncorrected n q . 

When experimental samples are not very large, fixed- iV correlation measurements become prob- 
lematic due to large sampling errors. Therefore one traditionally sums up all fixed- N cumulants 
121 Q, M to yield what are (inaccurately) called "short-range correlations" (SRC); in second 



order 

4 RC (xi, x 2 ) = {k 2 {x x , x 2 \N)} = J2 p N n 2 (x u x 2 \N) . (15) 

N 

Inclusive cumulants are then split up into the SRC contribution plus "long-range correlations" 
(LRC): 

k 2 {xi,x 2 ) = {k 2 (^i,x 2 |A^)} + {Ap N {x!)Ap N (x 2 )} , (16) 

where Ap N (x) = pi(x\N) — {pi(x\N)}. The SRC sum over fixed- iV cumulants is purported to 
represent an average over the standard fixed- iV correlations, while the second term represents 
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the LRC contributions resulting from the strong variation of p(x\N) with overall multiplicity N. 



As we have shown, it is preferable to replace the K q (x\N) by internal cumulants k^(x\N) to elim- 
inate the fixed- N multinomial contributions from the short range part. We therefore propose 
that the traditional SRC/LRC formula should be modified in favor of splitting the inclusive 
cumulants into iV-averages over "internal" and "external" correlations. In second order, this 
would be 

k 2 (x 1 ,x 2 ) = {k i 2 (x u x 2 \N)} (17) 
+ {^ u[t (x 1 ,x 2 \N)+Ap N (x 1 )Ap N (x 2 )}. 

LRC/SRC formulae up to fourth order were catalogued in Ref. [[?]]. In all cases, "internal" 
correlations would be given by {K q (x\N) — /t™ ult (a;|iV)}, while the "external" correlations cor- 
respond to the formulae in |7j plus the appropriate multinomial cumulant. 

For practical reasons, iV-averages of internal correlation measures over limited multiplicity 
ranges [A, B] are of importance, 



J2n=a Pn 



Fixed- N correlations are but one example of generalized cumulants. Equation (^) shows that 
subtraction of cumulants is appropriate when any kind of reference distribution is given. Ap- 
plications which come to mind immediately are subensembles characterized by a given number 
of jets and any Monte Carlo-generated simulation. Differences between real data and the dy- 
namics contained in a given Monte Carlo code would be quantifiable again by a generalized 
cumulant K = K data - k mc . 

Generalized moments p^ h are defined as functional derivatives of the ratio of generating func- 
tional, 

i ( _a \ /n,n \ 

; (19) 

A=0 



for the fixed- N case, the corresponding internal moments are related to the standard moments 
by 

p{(x\N) = 4(x\N) = Pl (x\N) , (20) 
pl{x u x 2 \N) = k i 2 (x u x 2 \N) + kKx^kKx^N) (21) 
= P2(x 1 ,x 2 \N) + (l/N)p 1 (x 1 \N)p 1 (x 2 \N) , 

and further, omitting the arguments for brevity, 

/ 3 2 \ 

p{ = P3+ Plplpl - M ' ( 22 ) 



K N N 2 , 

t 1 ( 5 6 \ 

Pa = P4+ff2^PiPiP*-PifrPiPi\Jf2-Jp)> ( 23 ) 

the brackets under the sum indicating the number of permutations. Internal moments have 
the property that when the measured moment p q (x\N) behaves multinomially, p q (x\N) — > 
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(N^/N q ) nf=i Pi(xi\N), they factorize without a prefactor: 

4(f|iV)^n Pl NiV). (24) 



For a subsample of fixed N, moments and cumulants at fixed multiplicity can be normalized in 
two ways. The most natural normalization procedure is to use in the denominator exactly that 
quantity which the numerator would default to if the process were fully independent. Hence 
the normalization for fixed- N factorial moments over some domain Q m would be 



F q (n m \N) 



) 

'.ZnJv\N"\ 

(25) 



In m dxpf^(x\N) 
N q J nm dxp q (x\N) 



ATM f nm dx Pl (x 1 \N)---p 1 (x q \N)' 

while the internal moment would be normalized according to 

t, ■ x fo dxpKxlN) , x 

f(i m dxpi(xi\N)---p 1 {x q \N) 

Both these definitions yield F q = = 1 for any Q m when the measured p q (x\N) behaves 
multinomially within the total window Qtot- For internal cumulants, the same normalization 
(p3) as for the internal moments would be appropriate. 



Internal moments integrate under the total phase space domain Q tot to J n dx p q {x) = N q , while 
p q (x\N) integrates to N^, so that both normalized moments become unity when integrated 
over Q tot for any distribution: 

F q I ({l tot \N) = F q (Q tot \N) = l. (27) 

The integral over Q tot of normalized internal cumulants will, of course, be zero just as the 
unnormalized ones. 

In summary, we have shown that generalized cumulants are an improved measure of correla- 
tions in samples of fixed multiplicity. The latter are more sensitive in discriminating between 
dynamical mechanisms and models than inclusive quantities. The procedure outlined here can 
be applied to a number of other conditioning problems within and beyond high energy physics. 
Detailed applications to correlation measurements such as pion interferometry are in progress. 
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